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L p BILINEAR QUASIMODE ESTIMATES 


ZIHUA GUO, XIAOLONG HAN, AND MELISSA TACY 


Abstract. In this paper, we investigate the L p bilinear quasimode estimates on compact 
Riemannian manifolds. We obtain results in the full range p > 2 on all n-dimensional 
manifolds with n > 2. This in particular implies the L p bilinear eigenfunction estimates. We 
further show that all of these estimates are sharp by constructing various quasimodes and 
eigenfunctions that saturate our estimates. 


1. Introduction 

Let (M, g ) be a smooth and compact Riemannian manifold without boundary. We denote 
A = A g the Laplace-Beltrami operator on M. An eigenfunction u of —A satisfies —A u = A 2 u 
with A its eigenfrequency. In 1988 Sogge [ Sol . So2j proved that 


Nb S 


< C\ 5< ' n,p ' > \\u\\ L 2 , where S(n,p) = 


Tn-rtl for2<p<2gnl, 

n —1 _ n 

2 p 


r 2(n+l) / / 

tor —— r 1 < p < oo. 

n —1 — 1 — 


(l.i) 


Here, C is independent of A. In 2007 Koch, Tataru, and Zworski [KTZj extended this result 
to quasimodes, i.e. approximate eigenfunctions in the sense that 

||(-A-A 2 H| l2 <CX\\u\\ L 2. 

In fact, their result holds for Laplace-like semi classical pseudodifferential operators. 

In this paper, we investigate bilinear eigenfunction estimates. That is, for two eigenfunc¬ 
tions u and v, we estimate ||mu||lp in terms of their eigenfrequencies. One can of course 


use Holder’s inequality and Sogge’s L p linear eigenfunction estimates (1.1) to prove a upper 


bound of ||uv||lp. For example, on a Riemannian surface (i.e. two dimensional Riemannian 
manifold), let u and v be two L 2 -normalized eigenfunctions with eigenvalues A 2 < /i 2 . Then 

IMU 2 < IMMIIMU 4 ~ As/xi - < p*, 

or with a different pair of Holder indices (among other possible choices) 

\\uv \\ L 2 < ||«||l°°||^||l 2 < At 

The second bound A5 does not depend on the higher frequency //, but is not necessarily 
better than the first bound (given, say, A ~ p). However in 2005, Burq, Gerard, and 
Tzvetkov [ BGT3] proved that 

IMU 2 < AT (1.2) 

This estimate is clearly better than both the two previous bounds. Moreover, they showed 
that A 4 bound is sharp on § 2 (see more discussion on the sharpness of bilinear eigenfunction 
estimates in Section [!]). 

This improvement is crucial in Burq, Gerard, and Tzvetkov’s investigation of nonlinear dis¬ 


persive equations on manifolds [ BGT3, IBGT4] , They used (1.2) to study the wcll-posedness 
of Cauchy problems involving nonlinear Schodinger equations on compact Riemannian sur¬ 
faces. In particular, they obtained the critical well-posedness regularity for cubic Schodinger 
equations on § 2 . 
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In fact, Burq, Gerard, and Tzvetkov [BGT3] proved (1.2) for spectral clusters. Subse¬ 
quently, they |BGT21 [BGT4j generalised (1.2) to higher dimensions as follows. 

Theorem 1.1 ( L 2 bilinear spectral cluster estimates). Let (M, g) be an n-dimensional com¬ 
pact manifold, and x £ Cq 
\l 2 = WqWl 2 = 1, we have 


IIxa(/)x m (^)II 


L 2 ~ 


X W~ A - 

A). Then for all 1 < A < /x and 

A 3 

if n = 2; 

A^ log A| ^ 

if n = 3; 

. n — 2 

A 2 

if n > 4. 


Moreover, all these estimates are sharp (modulo the log loss in dimension three). 

Note that if u is an eigenfunction with eigenfrequency A, then X\( u ) = u ■ So Theorem 0 
in particular applies to eigenfunctions. 

The bounds in Theorem 1.1 depend only on the lower eigenfrequency A. In this paper, 
we generalise Theorem 1.1 to L p bilinear eigenfunction estimates for p > 2. In this case, the 
bound depends on both of the eigenfrequencies. We are able to derive a full range of sharp 
(modulo some log loss when n — 3,p — 2) estimates in all dimensions. The two main themes 
that arise from our results are: 

(1). The sharp bound for the L p norm of the product of two eigenfunctions is better than 
that obtained by 


Holder’s inequality + Sogge’s L p eigenfunction estimates; 

(2). In the sharp bound, the higher eigenfrequency has a smaller exponent than the lower 
eigenfrequency. The extreme examples of this is L 2 bilinear estimates, the exponent of 
the higher frequency is 0. 

From a technical perspective it is just as easy to work with quasimodes (rather than exact 
eigenfunctions). So in this paper we prove the L p bilinear estimates for quasimodes. To this 
end, it is convenient to work in the semiclassical setting. Denote h = A -1 the semiclassical 
parameter. Then an eigenfunction u satisfies P(h)u(h ) = 0, where 

P(h) = -h 2 A - 1. 

We write p(x,hD ) as a semiclassical pseudodifferential operator with symbol p(x,f) and 
study 0 L i{h ) quasimodes as given by the following definition. 


Definition 1.2 (Quasimodes). A family {u(h)} (0 < h < h 0 1) is said to be an O^fih) 
quasimode of a semiclassical pseudodifferential operator p(x, hD ) if 

\p(x,hD)u(h)\\ L2 < h\\u(h)\\ L2 . 


In the semiclassical framework, we can treat quasimodes of semiclassical pseudodifferential 
operators that are similar to —h 2 A — 1 in the same fashion as quasimodes of the Laplacian. 
We call such operators Laplace-like (see Definition 2.2). Our main theorem states that 


Theorem 1.3 ( L p bilinear quasimode estimates). Assume that p(x,£) is a Laplace-like 
smooth symbol. Let u{h ) and v(o) be two families of Oj, 2 [h) and 0^2 (a) quasimodes of 
p(x , hD) and p(x , aD), respectively. Suppose that 0 < a < h < h 0 for /j 0 <C 1 and both u(h ) 
and v(cr) 


admit localisation property (see Definition 2.1) 
\\u(h)v(o)\\ LP < G nj ,(h,a 


Then 
u(h)\\ L2 \\v(a)\\ L2 , 


where for n = 2, 


{ , _i A_I 

h ±o 2 p 4 
h 2 p~ I a 2 ^~ i 

7-1 -- 1 

h 2(jp 2 


for 2 < p < 3, 
for 3 < p < 6, 
for 6 < p < oo; 
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3(n-l) I n+1 n —1 | n— 1 

h 3 a ~4 

n —1 1 | n 

h 2 cr 2 p 


2(n+l) 
n—1 5 


for 2 < p < 
for <p< 


oo; 


for n > 3, (n,p) 7 ^ (3,2), 

and, 

G 3t2 (h,a) = h~^\ log/?,|3. 

Moreover, all these estimates are sharp (modulo the log /oss m t/ie case (n,p) = (3,2)). 

is an O^(h) quasimode of —h 2 A — 1 


1.1 


We point out that Xx(f) con sidered in Theorem 
with h = A -1 . (See Section 5.4 ) Similarly, x^(g) is an 0 L 2(a ) quasimode of — a 2 A — 1 with 


cr — p 1 . Therefore, the L p bilinear spectral cluster estimates follow as a simple consequence 
of the above theorem. Moreover, they are also sharp. 

Theorem 1.4 ( L p bilinear spectral cluster estimates). Suppose that p > X > 1. Then 

\\X\(f)x^9)W L 2 < G n , p (\-\p- u 


L 2 l|y||L 2 , 


where G n ^ p is as given in Theorem 1.3. In particular, let u and v be two L 2 -normalized 


eigenfunctions with eigenvalues X 2 and p 2 . Then 


\\uv\\ LP G n , p (A , P )• 

Moreover, the estimates are sharp on the sphere. 

Application of bilinear eigenfunctions to nonlinear dispersive equations. An im¬ 
portant application of the bilinear eigenfunctions estimates lies in the study of nonlinear dis¬ 
persive equations on Riemannian manifolds. Consider the Cauchy problem involving cubic 
(focusing or defocusing) nonlinear Schrodinger equation on a compact Riemannian manifold 


id t u + An = ±|w| 




(1.3) 


w(0, x ) = wo(t). 

There is a long history in studying the well/ill-posedness of this system, see e.g. 

Gerard, and Tzvetkov [ BGT2 , BGT4 J and the references therein. Denote II 
space on M. Burq, Gerard, and Tzvetkov [ BGT2 , Theorem 1] proved that (1.3) is locally 
well-posed in H s for s > 1/4 on the two-sphere § 2 . Their argument combines the L 2 bilinear 


Burq, 
the Sobolev 


eigenfunction estimate ( 1.2[ ) with the eigenvalue distribution on § 2 . On the questions of 
ill-posedness, their earlier work [ BGTl j showed that (1.3) is ill-posed in H s for s < 1/4. 
Therefore, H 1 ^ is the critical regularity on § 2 . We express this by saying that the critical 
threshold s c (§ 2 ) = 1/4, following the notation in [BGT 2 |J. 

It is only on special manifolds (such as the sphere) that such critical regularity threshold is 
known. Another known example is the torus. Bourgain [B] and Burq, Gerard, and Tzvetkov 
jBGTIj . showed that s c (T 2 ) = 0 where T 2 is the two-dimensional torus. It remains open to 
find the critical threshold on general Riemannian manifolds, and to relate it to geometry of 
the manifold. For more information about nonlinear Schrodinger equations on Riemannian 
manifolds, we refer to [B, BGTl , BGT3 BGT4 J and the references therein. 


Related literature. Burq, Gerard, and Tzvetkov [ BGT4 , Theorem 3] also proved the fol¬ 
lowing If 2 trilinear eigenfunction estimates in dimension two and three. 


Theorem 1.5 (L 2 trilinear spectral cluster estimates). Let (M, g) be an n-dimensional com¬ 
pact manifold and x £ Co°(M). Given 1 < A < p < u and \ |e11x , 2 = ||/||l 2 = Ilfi'IU 2 = 1, we 
have 

ifn = 2 ; 
if n = 3. 

Here, £ G (0,1]; C £ depends on e and M; but C and C e are independent of X, p, and v. 


\\X\( e )Xn(f)Xv(g)\\ l 2 ^ { ^\i^s L—e 

p 2 
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Remark. We point out that when // « 
derived from the above L 2 trilincar one. 


A our L 3 bilinear estimate in Theorem 1.3 can be 


See also the recent work of Yang [YJ on L 2 multilinear eigenfunction estimates and their 
applications to a more general class of nonlinear Schrodinger equations. 

We remark that the bilinear eigenfunction estimates are sharp on the sphere. However, 
they are far from optimal in the case of torus. For example, on the two dimensional torus 
T 2 , given two eigenfunction u and v with eigenvalues A 2 and /i 2 , 

\\uv \\ L 2 < ||w||z,4||i;|| i 4 < C, 

where C is independent of A and //. This follows the classical result of Zygmund |Zy| that 

IMU* <c. 

Organisation. We organise the paper as follows. In Section [2j we use semiclassical analy¬ 
sis and localisation to reduce the problem to a local one. In Section [3j we prove Theorem 


priate quasimodes and spherical harmonics. Section [5] is an appendix containing the basic 
semiclassical analysis on which this paper relies. 

Throughout this paper, A < B (A > B) means A < cB [A > cB ) for some constant c 
independent of A or h; A za B means A < B and B < A] the constants c and C may vary 
from line to line. 


1.3 In Section Hi we show that the L p bilinear estimates are sharp by constructing appro¬ 


2. Semiclassical analysis and localisation 


In this section, we use semiclassical analysis to reduce the L p bilinear quasimode estimates 
to a localised problem. For the readers’ convenience, we include a brief discussion on the 
background of semiclassical analysis in the appendix (Section [5]). 

As any two quantisation procedures differ only by an O^lfi) term, we adopt the left 
quantisation for symbols, that is, 


p(x, hD)u(x) 


1 

(27rh) n 


e* <x v,Z)p(x,£)u(y)d£dy. 


We require that both functions u(h) and v(a) are semiclassically localised with respect to the 
relevant parameter h or er, respectively (see Definition 2.1) and that p(x,£) is Laplace-like 
(see Definition |2.2 ). 


Definition 2.1 (Semiclassical localisation). A tempered family {u(h)} is semiclassically lo¬ 
calised if there exists a function y G C'o°(T*M) such that 

u(h) = y(x, hD)u(h ) + Os{h°°). 


An immediate consequence is that if a family u = u(h ) is localised then 

\\u\\ Lq < Ch n (v~ 


- n[ ~v)\\u\\ LP + 0(h c 


( 2 . 1 ) 


where 1 < p < q < oo (see Theorem 5.5). In particular, 


_ n | n 

M LP <h 2 p \\u\ 

L -i\ 


L 2 • 


Note that the eigenfunctions Uj = u (A ■ ) admit localisation property. However, the above 


inequality gives a worse estimate than Sogge’s estimates in (1.1). 


Definition 2.2 (Laplace-like symbols and operators). We say that a symbol p(x,£) (or its 
associated semiclassical pseudodifferential operator) is Laplace-like if 

(1) for all (x 0 ,£o) e T*M such that p(x 0 ,£o) = 0, V^p(x 0 ,^o) 7^ 0; 

(2) for all x 0 £ M, the hypersurface {£ G T *M : p(x 0 ,£) = 0} has positive definite second 
fundamental form. 
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Now suppose that u(h) and v(a) are L 2 normalised semiclassically localised quasimodes of 
p(x,hD ) and p(x, crD), respectively: 

p(x, hD)u(h) = 0 L 2(h), and p(x,aD)v(a) = 0 L 2(a), 

where 0 < a < h < ho. From the assumption, they both admit the localisation property. We 
may assume that 

u(h) = x ( x ) hD)u(h) + Os{h°°), and v(cr) = x( x , crD)v(a) + O^u 00 ), 
where \ ^ C%°(K) with K d T* M. Hence, 

\\uv\\lp = \\x( x ,hD)u(h)x( x ,crD)v(a)\\ LP + 0(h°°)||x(a;, aD)v(a)\\ LP 
+0(a°°) ||x(x, hD)u(h) || LP + OihTa 00 ) 

< \\x( x ,hD)u(h) X ( x ,crD)v(a)\\ LP + 0(h°°)a- 5{n ' p) + 0(a°°)h- 5{n ’ p) + Oih^a 00 ) 

< \\x{ x ,hD)u(h)x( x ,crD)v(a)\\ LP + h N a~ 5M + a N h~ 5 ^ p \ 

Inspecting the definition of G ntP (h,c r) and 5(n,p), we see that a~ 5(jl,p ’ < G niP (h,a). Hence, 
we only need to estimate the first term. Now we further localise each quasimode. Let 

N 

X = ^Xi, 


i— 1 


where N < oo and Xi l ias arbitrarily small support. Thus, 

N 


\\x( x ,hD)u(h)x(x,aD)v(cr)\\ LP = 


N 


Xi( x , hD)u(h) Y Xj( x , crD)v(cr) 


i =1 
N 


3 =1 


LP 


< Y || Xi( x ,hD)u(h) ■ Xj( x ,<?D)v(<r)\\ LP . 
i,j =1 

We then reduce the problem to a local one by separately estimating each term 

\\Xi(x,hD)u(h) ■ Xj( x ,vD)v(a)\\ LP 

in the summation. Since semiclassical pseudodifferential operators commute up to an 0(h) 
error (see Section |T3 ), we have 

p(x, hD)xi(x , hD)u(h) = Xi( x , hD)p(x, h,D)u(h ) + hr(x, hD)u(h). 

Therefore, if u(h) is an 0 L 2(h) quasimode of p(x, hD ), then Xi( x i hD)u(h) is also an 0 L 2(h ) 
quasimode. 

In each patch, we may assume the Riemannian volume coincides with the Lebesgue density 
in local coordinates. Since the localisation property is consistent with 0 12 (h) quasimodes, 
we have 


p(x, hD)xi(x, hD)u(h) = 0 L 2(h), and p(x, aD)xj(x, aD)u(a) — Ol 2 (<t)■ 

We may further reduce to studying a product Xi( x i hD)u(h)xj( x , crD)v(a) where p(x,£) is 
zero somewhere on supp XjUsupp Xj- Suppose this were not the case, then one of the following 
two cases happens. 


2.1. |p(x,£)| > c on suppxj. Since 

p( x i hD)xi( x , hD)u = hf, where ||/| l2 < ||u|j i2 , 
and p(x, hD) is invertible on the support of Xh we have 

Xi( x , hD)u = hp(x , hD)~ l f , 
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This implies by (2.1) that 

\\Xi(x,hD)u\\ LOO < h-h~% \\u \\ l2 

because Xi( x ,hD)u is semiclassically localised. Therefore, applying Holder’s inequality and 
Sogge’s L p estimates fll.l ) on Xj{ x , crD)v, we have 


\\Xi(x,hD)uXj(x,aD)v\\ LP < \\xi(x,hD)u\\ LOO \\xj(x,aD)v\ 

< h-f+V- 5 ^) ||n || L 2 ||n|| i2 , 


Lp 


which is a better estimate than given by Theorem |1.3 
2.2. |p(:r,£)| > c on suppxj. Then 

\\Xj(x,aD)v\ L2 <a\\v\\ L2 . 


This implies by (2.1) that 


\\Xj(x,aD)v\\ Lr <a-a n 0 2 ) |n| i2 


because Xj( x ,crD)v is semiclassically localised. Therefore, applying Holder’s inequality and 
Sogge’s L p estimates 0 on Xi( x , hD)u , we have 

I Xi (x , hD ) uxj (x , aD) v || LP < \ \ Xi (x , hD ) u \ \ L oo \ \ Xj (x, crD) v \ \ LP 

< h 2 a h h H| L2 |M| L 2, 

which is also better than that of Theorem 11.31 

Therefore, we assume that there is an (xo>£o) e su PPXi and an (5)o,£o) G supply such 
that 

p(x o,£o) = 0, and p(x o,£o) = 0. 


3. Proof of L p bilinear quasimode estimates 
Following the reduction in the previous section, we may now work with the product 

Xi(x,hD)u-Xj(x,aD)v, 

where x.% and Xj have small support, and there is an (£ 0 ,£o) £ suppxi and an (iio^o) £ 
supp Xj such that 

p(x 0 ,£o) = 0, and p(x 0 ,£ 0 ) = 0. 

In these coordinate patches, we often drop Xi and Xj and for notational convenience simply 
write u and v. At this point we use our assumptions on p(x , £) to factorise the symbol. 
Symbol factorisation will allow us to write 

P(x,0 = e(x, £)(£* - a(x,£')), |e(x,£)| > c 

for some £j. We associate Xi with time t and treat u as an approximate solution to an evo¬ 
lution problem. We then obtain L p estimates by considering them as Strichartz estimates 
with equal weighting on space and time. This “variable freezing” is a well-established tech¬ 
nique; see for example Stein [St j , Sogge [ Sol . 15o2|. Mockenhaupt-Seeger-Sogge [ MSB ], and 
references therein. In the semiclassical setting Koch-Tataru-Zworski [ KTZ J used it to obtain 
L p estimates for quasimode in a similar setting to that we are considering. 

Recall that the first assumption for p(x, hD) to be Laplace-like is that whenever p{x 0 , £ 0 ) = 
0 the gradient V^p(x o,£o) 7 ^ 0. Therefore, there is some £fc for which d^ k p(x o,£o) 7 ^ 0. In 
fact, we may choose a £& (which we will call £ 1 ) so that if £ = (£i,£'), then 

V € 'p(x 0 ,£ 0) = 0. 

We now divide our problem into two cases. 

Case 1. 


d^p(x 0 ,£o) 0; 
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Case 2. 

d^p(x o,£o) = 0. 

Case 1: By the implicit function theorem, we may factorise the symbol p(x,£) such that 

p(x,0 = ei(x,0(6 -ai(a;,0) on su PPXu 
P(z,0 = e 2 (x,0(£i~ a 2(x,O) on suppxj, 

where both \e 3 {x, £)|, \e 2 (x, £)| > c > 0. Hence, we may invert e\(x,hD) and e 2 (x,aD) to 
obtain 

(hD xl -ai (x,hD x '))u(h) = 0 L 2(h\\u(h)\\ L 2), 

and 

(crD Xl - a 2 (x,aD x >))v(a) = Oi?{a\\v{a)\\v). 

We will treat this case in Propositions |3.1| and |3.2| by setting X\ — t and reducing the problem 
to one involving semiclassical evolution equations. Picking np on that notation we say that 
in Case 1, u and v propagate in the same direction. Note that 

{£ G T*M : p(x o,0 = 0} = {£ G T * x M : & = ai(x 0 ,O}, 

and 

{£ e T*M : p(x o,0 = 0} = {£ G T*M : 6 = a 2 (xo,0}- 

Now since 

V^ai(x 0 , £o) = 0, 

the second fundamental form is given at this point by 

d 2 a\ 

so since p(x, £) is Laplace-like the matrix, 

d 2 d 1 

is positive definite. A similar argument gives that 

d 2 a 2 

is also positive definite. 

Case 2: Since d^p(x 0 , £o) = 0 but V^p(x 0 , £ 0 ) 7^ 0, there is some other (which we will call 
£ 2 ) such that 

<%pOH),£o) ^ 0. 

We choose this direction so that if £ = (£i,£ 2 ,£), then 

V(p(x, 0 ,io) = 0. 

Now we write 

p(x,£) = ei(x,£)(fi - ai(®,6,0) 011 su PPXu 
p(x,0 = e 2 (x,£)(£ 2 -a 2 (a:,£i,£)) on supply. 

Again, |ei(x,£)|, |e 2 (x,£)| > c > 0. So we invert e\(x,hD) and e 2 (x,aD) to obtain 

(hD Xl - ai(x, hD X2 , hD x ))u = 0 L 2(h\\u(h)\\ L 2), 

and 

(aD X2 -a 2 (x,crD xl ,aD x ))v = 0 L 2(a\\v(a)\\ L 2). 
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We will treat this case in Proposition |3.3| by setting x\ = t\ and x 2 = t 2 and reducing the 
problem to one involving semiclassical evolution equations propagating in different direc¬ 
tions. We adopt the notation x = (xi,x 2 ,x) and £ = (£i,£ 2 ,£). Note that by choosing the 
support small enough we may assume that 

|(<% 2 ai, Vfdi)| < e, and |(%a 2 , V f a 2 )| < e. 

As in Case 1, we obtain that the matrices of second derivatives <9^ and <9^ ^a 2 are 
non-degenerate. 

We are now in a position to complete the estimates of Theorem 1.3| by studying the above 
two cases. 

3.1. Propagating in the same direction (Case 1). We split this into two further cases. 
Case la when a < ch (for some small constant c dependent only on the manifold), in this 
case the frequencies of the quasimodes are significantly different. Case lb is of course when 
ch < a < h, in this case the frequencies are close enough to be treated as the same. 

Proposition 3.1 (Case la). Assume the notations in Case 1: 

IK hD Xl — afix, hD x i))u(h )\\ L 2 < h\\u(h)\\ L2 , 

and 

II0 D X1 - a 2 {x, aD x ,))v(a)\\ l2 < a K<x) || L2 , 

where d|,^,ai and d|,^,a 2 are positive definite. Suppose that a < ch for some small c dependent 
only on the manifold. Then 

\\uv\\lp ~ G n , P (h, cr) \\ u \\l 2 Ill’ll L 2 ’ 
where Gn^h^a) is given in Theorem 1.3. 

Proof. We adopt the notation t = x i and write x G M as x = (t, x'). Let 

E h [u] = ( hD Xl - ai(t,x', hD x t))u(h), 

and 

E a [v\ = ( aD xl - ai(t,x / ,aZ9 ;c /))n((j). 

Using Duhammel’s principle, we have 


and 


where 


and 


u{t,x') = U h (t, 0)w(0,a/) + - / U h (t - s,s)E h [u\ds, 


v(t, x') = U a (t, 0)u(0, x') -\— / U a (t — s, s)E a [v] ds, 

v J o 

(hD t - afis + t, x ', h.D x f))U h (t , s) = 0, 

U h (0,s) = Id, 

(aD t - a 2 (s + t, x\ aD x ,)U a (t, s ) = 0, 

U a (0, s) = Id. 


For simplicity of notation, we write £4(f) := Uh(t,0),U a (t) := U a (t, 0). We may then write 
uv as 

1 r f ~ 

u(t , x')v(t, x') = Uh{t)u(0, x / )U cr (t)rK0, x') -|— Uh(t)u(0, x') / U a (t — s, s)E a [v] ds+ 

(T .In 


1 


t rt 


h 


ha 


-U a (t)v( 0,x') / U h (t-s,s)E h [u\ds+— / / U h {t-s 1 ,si)E h [u\U a (t-s 2 ,S 2 )E a [v\dsids 2 . 


0 J 0 
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Since t is localised and 

\\K[v]\\<cr, \\E h [v]\\ L2 < h, 

it suffices to prove that 

U h (t)f(x')U a (t)g(x' 


< 




Gn,p(h,°) \\f\y, hli?, • 


(3.1) 


Using the parametrix construction in Section 5.5, we have 




and 

where 

and 


U a (t)g(x') — ——[ I e- (02(ty ’ ,/) {w '’ 1,,)) b 2 (t,x',ri')g(w')dw'dr )', 
d t <f>i(t,x?,£') - a 1 (t,ar', V x /0i) = 0, 0i(O, x 1 , £') = (x',£'), 


(3.2) 

<%0 2 (f,a/,7/) - a 2 (t,x', Vx'fc) = 0, 0 2 (O,x',t/) = (x',r)'). (3.3) 

Therefore, we need to study the L 2 (M 2n ~ 2 ) —y L p (M. n ) mapping properties of the operator 
W Kc (t) given by 

W ht<r (t) = ^ n-1 J 6 ( t> dTy'de'ds/dw', 

where 

K^'^W) = b 1 {t,x , ,Z')b 2 (t,rf,r/). 

We will view this as a Strichartz type estimate with equal weighting on space and time. For 
unitary operators V(t), such estimates depend only on the dispersive bounds 

lV(t)V*(a)l Ll ^ 00 <|t* a |-". 

See for example the abstract framework given in Keel and Tao EH. Our operator Wh, a (t) is 
not unitary but we may still use this framework as in Tacy m by proving L 2 —* L 2 estimates 
to replace unitarity. We prove 

|^(«>)IL. -+l« S + 1 1~ s|)^»(<t + 1 1- 

X 1 X 1 

and 

WWkAWiMWv < h-”a-»*{h + 1 1- s|) _1I 0 + \t - s\)-«. 

x f x f 

Then by interpolating between then obtain 

l|WV(t)WS>)||y-_y $ h-°’a-e’(h + |t - s|)^(a + \t - a|)-">. 

x r X 1 

As in the Keel-Tau framework and in earlier work by Mockenhaupt, Seeger, and Sogge [MSS] 
and Sogge isn. we use Young and Hardy-Littlewood-Sobolev inequalities to resolve the 
1 1 — s | integral. Therefore, we need to calculate Wh, a (t)W^ a (s), where 


Here, 

W ( t , x\ s, z') = 


W hia (t)W^ a {s)f(x') — / W(t, x', s, z')f (z')dz 1 . 


a> ^ v ' llV f 2 ) dA 

dA = dy' dw' d^d^drj'idrjz, (3.4) 


(/ra) 2 0 d 


where 


^h,a(t,s,x ,y ,w ,z ,£i) “ M S > Z ’Q ~ (y ,€i ~ £2) 



















10 


ZIHUA GUO, XIAOLONG HAN, AND MELISSA TACY 


a 


02 (t, x', T) i) - 01 (s, z'l r)’ 2 ) - (y 1 , T)[ - T] 2 ) 


We use stationary phase asymptotics [ Zw , Theorem 3.16] to calculate the (]/, £[, w', rf -]) inte¬ 
gral in (3.4). As the stationary point is always non-degenerate, we obtain 


x-- / s, a/, 2 / ?/) drf 

a n 1 j 

:= K h (t,x',s,z') ■ K a (t,x',s,z'). 

We analyse K h and K a by considering critical points of 

0i(f,®',O -0i(s, *',£') in 

and 

x\ rf) - 0 2 (s, z\ rf) in 77 ' 

To that end, following [KTZ j. we study the oscillatory integral 




(3.5) 


where 


V>i(M, s', *',£') = 0i(t,x',^') - 0 1 (s,z',f')- 

The critical point of 0i solves 

(t,s,x',z',^) = V e [ 0 i(t,x',^) - 0 !(s,^,O 


= 0. 


(3.6) 


From (3.2), we have that 

0i(f,z',O - 0i( s > d, £) = (x' - + 0(|s|)) + (t - s)(ai(s, 2 ; / ,^) + 0(\t - s|)). 

Then (3.6) implies that 

0 = (1 + 0(|s|)) ■ (x 1 - z) + (t — s)(d^a 1 + 0(\t - s|)), 
and thus for critical points to occur we must have 

\x' — z'\ = 0(\t — s|). 


If \x 1 — z'\ \t — s|, using integration by parts to (3.5) we get for any iVeN 

I K h {t,x\s,z')\ < ( 1 + ^“^) • 

If \x' — z'\ < |t — s|, there might be critical points. The Hessian is given by 

<9|,£,0i(t, s,x',z' : £') ={t - s)(d^,ai + 0(\t - s|)) + 0(\x' - z'\s) 

=(t — s)(d^,ai + 0(\t\ + | s |)). 

Thus by the non-degenerate condition there exists at most one critical point. We know that 
<9|,£,a! is a non-degenerate matrix, so 

|det (9|, ? /0i) | > c\t — s| n_1 . (3.7) 


A similar calculation holds for 


where 


K a (t, x', s, z') = — I e^ t ’ s ’ x '’ z '^)b 2 (t,s,x',z , , V ')dri', 
0 2 (t, s, x', z', rf) = 0 2 (t, X 1 , rf) - 0 2 (s, z', rf). 


(3.8) 
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If \x' — z'\ \t — s|, then 


\K a (t,x',s,z')\ < 


n— 1 


a 


1 + 


\x — z 

<J 


-N 


If \x' 


< 


\t — s|, there exists at most one critical point that satisfies 


V r? /'0 2 (^, s, x\ z\ rj') 




Mtix',r)') - (j) 2 (s,z\r)') 


The Hessian satishes 

|det ip 2 ) \> c\t - s| n_1 
We split the kernel into three parts 

• Regime 1. \t — s| < Ca, 

• Regime 2. Ca < \t — s| < Ch, 

• Regime 3. Ch < \t — s|, 

for some C independent of t, s, h, a. We write 


0 . 


(3.9) 

(3.10) 


W(t, xs, z') = W x {t, x', s, z') + W 2 (t, x', s, z') + W 3 (t, x', s, z'), 
where x e Cg°([-1,1]) and 

Wi(t,x',s,z') = x((T _1 |t - s\)W(t,x',s,z'), 

W 2 (t,x',s,z') = (1 - x(cr _1 |t - s\))x(h~ L \t - s\)W(t,x',s,z '), 

W 3 (t,x',s,z') = (1 - x( a ~ 1 \t - s\) - (l-xi^lt ~ •s|))x(/i“ 1 |t - s\))W(t,x',s,z'). 

In the following computation, we slightly abuse the notation and use Wi(t, s ) as the operator 
with integral kernel W t (t, x', s, z'), i = 1, 2, 3. 


3.1.1. Regime 1. We have that \t — s| < Ca in the support of W\. In this case, note that as 
there cannot be a critical point in either £' or rj if \x' — z'\ > C\t — s[, we can integrate by 
parts in rj to obtain 

|Wi(t,z , ,s J z , )l < h-^-Va-^-V (\ + 



for any natural number N. Therefore, applying Young’s inequality (see e.g. [So3 . Corol¬ 
lary 2.1.2]) we obtain that if W\(t, s) is the operator associated with the integral kernel 
W\ (t, x', s, z'), then 




2(n — 1) 
V 


(3.11) 


3.1.2. Regime 2. We have that Ca < |t — s\ < Ch in the support of W 2 . If \x' — z'\ \t — s\ 
or K a has no critical point, then we immediately get the same bound (3.11) and hence we are 
done. Now we assume K a has a unique critical point rj c . By the stationary phase theorem 
m Theorem 7.7.5], using (|3.8|) and (|3.10[), we have 


W 2 (t, s,x ', z') = Kh(t, s,x ', z') ■ a 2 \t — s| 2 e° e ^' s ' x ’ z ^B(t, s, x', z'), 


where 


Q(t, s, x 1 , z') = i/j 2 (t, s, x', z', rj c ) = <j) 2 (t, x', g' c ) - <j> 2 (s, z', 77 '). 
Thus we immediately get 

\\W 2 (t,s)\ Ll L oo < h-^-^a-^lt - s\ 

x' x' 

To get the L 2 , —> L 2 , estimates we write 


re —1 
2 


(3.12) 


W 2 (t,s,x',z') = / W 2 (t,s,x',z',£')d£', 
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where 

W 2 (t, s, x', z ', £') = ^ f )#( t> s> x ' ; -'). 

Due to the compact support in £' it is enough to estimate the L 2 —> L 2 norm of W 2 (t, s,£'), 
the operator with kernel Wiit, s, x', z', £'). 

We will do this by writing 


||W 2 (M,0/llk = I W 2 (t,s,y',z',e)f(z’)f(y')dy’dz' : 


IIL J , 


where 


W 2 (t,s,y', = / W 2 (t,x',s,y',(')W 2 (t,x',s,z',£')dx'. 


(3.13) 


and estimating |hW 2 (t, s,y', z')\ by integrating by parts in x'. To do this we need to know how 
B(t,s,x',z ') behaves under differentiation in x'. 

If we differentiate (3.9) in x’ gives 


0 = Id + 0(\t\ + \s\) + (t- s) ([d x ^ c ] T d^^cj ) 2 + 0(\t - s|)) . 


Hence, 
and thus, 


-1 


d x >r)' c {t,s,x',z') = - -(Id + 0(\t\ + \s\)){d 2 , ,fa) \ 


t — s 


(3.14) 


\d Xi ri' c (t,s,x',z')\ < 


C 


We can repeat the process to get 

I Dl,r)'c{t,s,x',z')\ < 


1 1 — s\ 


C 

1 1 — s\ 


M 


for any multi-index 7 . Furthermore, we know from Tacy m Lemma 4.1] that 


\Dl,B(t,s,x',z')\ < \Dlr]' c (t, s, x', z')\ < 


C 


1 1 — s\ 


It! 


(3.15) 


We may now proceed in integrating (3.13) by parts in x’. 

The phase function is 

da,hit, s, x', y', z') = ^ [9(t, s , x', y') - 0 (t, s, x', z’)\ + ^ [ 0 i(s, y\ £') - fa (s, z', £')]. 


Note that the term 


^ [fa{s,y',0 ~ fa(s,z',£')]. 


does not depend on x'. Expanding with a Taylor series about y' = z’ we obtain 
8(t, s, x 1 , y') - 6(t, s, x\ z!) = ^ V z '0(t, s, x\ z') ■ (y - z!) + Ofay - z'\ 2 ) 
So 

^ x ’K,h{t,s,x\y',z') = -\d 2 x , z ,9{t,s,x\z') ■ {y' - z') + 0{\y' - z '| 2 ) 

a L 

Therefore, we need an expression for d 2 , y ,9. We have that 

0(t, s, x\ z') = fa(t, x\ 77 ') - fa(t, z\ 77 '), 

and 

0 = d v '{fa{t,x',r/ c ) -fa(t,x?,r/ c )). 

Differentiating (3.16) in z' and x' we obtain 

d 2 e 


d 2 fa drf c drf c 
dxidzj dr)idr) k dx t dzj ’ 


E 


(3.16) 
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which in matrix form is 


d 2 e 

dx'dz' 


= (9.A] 


(3.17) 


We already have an expression for d x ir]' c (t, s, x', z') from (3.14). A similar process gives us 
d Z 'rf c (t,s,x',J) = ——— (Id + 0(\t\ + |s|))(^ v 0 2 ) _1 . 

L S 


Combining this with (3.17), we have 

o 2 e 


dx'dz' t — s 


(9,v^2 +0(|«-s|)). 


(3.18) 


Therefore, 


|V x #0 tr (t,s,a; , ,2/ , ,z , )l > 


c\y' ~ z'\ 
<j\t — s| 


So each integration by parts of (3.13) with respect to x' gains a factor of 

a\t — s| 

W ~ A' 


but, in view of (3.15), looses a factor of \t — s| 1 from hitting the symbol. Overall each 

a 


integration by parts gains a factor of 


Therefore, we obtain 

\W 2 (t,s,y',z',^)\ 

< h-^ n -d a -^-d 1 1 - (l + 


W - A 


W - A 


a 


-N,. 


X 


W - A 

1 1 — s\ 


X 


| z — x 
1 1 — s\ 


dx' 


</r 2( " _1) <7- (n-1) ^1 + ^ ^ ^ 


— N 


So by Holder and Young inequalities, 

iw2(i, s .r)/ih<ft- 2( ”- 1) i/ii,. 

X' x' 

which implies 

x' x' 

and 


„ < L-O-l) 

r 2 _. r 2 > 


|W 2 (i,s)|| 


x' x' 




(3.19) 


3.1.3. Regime 3. We have that \t — s\ > Ch in the support of W 3 . If \x' — z'\ \t — s\ or K a 
has no critical point, then W 3 has the same bound as W\ given by (3.11) and we are done. 
Now we assume \x' — z'\ < \t — s| and K a has a critical point r]' c . If Kf, has no critical point, 
thenlYi has the same bounds as W 2 given by (3.12) and (3.19) . So we may also assume K h 
has a critical po: 

W 3 (t,x',s,z') = 


also has a critical point, denoted as £'. Therefore, we obtain 

1 


SC 

— | t- s\-^ 11 - s\-^ e K^Rt,s, x ',z')+^ 2 (t,sy,z')) B ^ a y j z '^ 


(ha) 


where 

and 


ei(t,s,x',A = ~(j)i(t,x',Q, 

e 2 (t, s, x\ z') = <j) 2 (t , x', r\ c ) - (j) 2 (t, x\ r \ c ) 


for £' and rf c the solutions to (3.6) and (3.9), respectively. Further we still have 

I D l' z 'B{tiSix',z')\ < 1 1- s| H71 
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for any multi-index 7 thanks to (3.15). From this representation, we can directly read off the 
L x , —» L™ estimates, that is, 


\\W 3 (t,s)\\ Ll ^ <h n z a ' l 2 \t-s\ (n 1} . (3.20) 

x' x' 

We need only then calculate the L x , —>■ L x , norm. As in Regime 2, we have 

112 


\\ W 3 (t,s)f\\ L 2 i = / W 3 (t,s,y',z')f(y')f(z')dy'dz'. 


Here, 

W(t,s,z',y) 


in which 
and 


=j W(t,x',s,y')W(t,x',s,z') dx' 

=/r- (n - 1) a- (?l - 1) |f - s\- 2{n ~ l) J s> x \ z ')B{t, s, x', y') dx', 

0 i(t, s,x',y',z') = 0 1 {t,s,x',y') - 6 V (t,s,x',z'), 


6 2 {t, s, x', y', z') = d 2 (t, s, x ', y') - d 2 (t, s, x', z'). 

Note that as cr is much smaller than h the major oscillatory terms come from this parameter. 
We need to calculate 

V x /(0i(t, s,x',y',z') + 0 2 (t,s,x' ,y' ,z')) 

Again we expand around the point y' = z! and obtain 

^x’( 0 i(t,s,x',y',z') + 6 2 (t,s,x',y',z')) = (^dl, z ,Qi + ^d 2 x , z , 6 ^\ • {y' - z') + 0 (\y' - z !| 2 ). 


We already know from (3.18) that 

d 2 x -z' e 1 = ^(^r 1 + 0 (\t - s|)), and d 2 x , z , 6 2 = -^-(d 2 ,^ 1 + 0 (\t - s|)). 
Since h > Ca, we have 


VAO 1 + O 2 ) 


> 


c\y' — z'\ 


a\t — s | 

So we may proceed in the same fashion as in the second regime to obtain 


I y' - A 


cr 


-N 


X 


Iiw , 3(«, s )/||k<r("- 1 >|i- s r(”- 1 >||/|l7, 


I w 3 (t, a,y\z ')\< |«- s r 2 <"-» ( 1 + 

This implies that 
and 

|W' 3 (m)| j 

Recall that we needed to prove 

\\(W h At)q)(x')\\ LP < G(/r,a)||g(i/ , ,w , )|| L 2, ,. 

t,x' y 

By TT* argument it suffices to show 


X 


\z' — x'\ 


\y' - A 

\t — s| J V V 1 1 — s 


dx'. 


. 7 _ re—1 1 | _ n—1 

II ~ h 2 I* - S l 2 • 

X' X' 


(3.21) 


Wi{t,s)f(s) ds 


< G(h , af 


L i*' 


L L> 


i = 1,2,3. 
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Regime 1: \t — s\ < Ca. We get from (3.11) and Holder’s inequality that 


Wi(t, s)f(s) ds 


< 


LP t y 


h (n 1] (T (n 1)+ 1 \\f(s, -)|U' ds 


\t—s\<C(j 


L p , 


L p t 


2 n 
P 


II Lf , ' 


Regime 2: Co < \t — s| < Ch. Interpolation between (3.12) and (3.19) gives that 

lW 2 (t,s)\ L ^ L , 

X f X' 

Hence, from Holder’s inequality, we have that 


/ i \ n —1 i n — 1 , , n — 1 . n — 1 


W 2 (t,s)f(s)ds 


< 


L L' 


t -i \ n — 1 i n — 1 , , n—i . n—i .. 

hr^-^o 2 + P \t- s \ 2 P | f(s,-)\\ p ’ds 


n—1 | n—1 




L p , 




2 n 
P 


L L' 


when 


P ~f~ 


2(n + l) 
n—1 


At p = ’yyvp, we use Hardy-Littlewood-Sobolev inequality in dimension one to get 


W 2 (t,s)f(s)ds 


< 


q x , 


/ n — 1 i n—1 


n — 1 | n — 1 


a ~\t— 


\t~s\ 2 p ||/(s,-)| J / # d« 


q 


/ \ n — 1 , n — 1 

< 3-+-ir 


" Lf , ' 

t,x' 


Regime 3: \t — s| > Ch. Interpolation between (3.20) and (3.21) gives that 

X f X 1 

Hence, from Holder’s inequality, we have that 


W 3 (t,s)f(s)ds 


< 


L L' 


, n — 1 n—1 , n—1 i (n El n ~ 1 

h~^2~a 2 ' r> \+ — o\ \ n r> 




3(n—1) I n+1 _ n — 1 ■ n — 1 

<h 2 p cr 2 p 


1 1 — s | 


L ty 


p \\f(sr)\\ LP 'ds 




when p (n + l)/(n — 1). But (n + l)/(n — 1) <2 unless n = 2,3. So we only need to 
consider the cases (n,p) = (2,3) and (n,p) = (3,2). 

For (n,p) = (2,3), by Hardy-Littlewood-Sobolev inequality in dimension one we get 


W 3 (t,s)f(s)ds 


< 


L \ x , 


\t~ 


h ^o ®|t-s| 2/3 \\f(s, -)l L 3/2 ds 


q 


<h 2(7 S ||/|| l 3/2 
t,x f 


For (n,p) = (3, 2), we can not use Hardy-Littlewood-Sobolev inequality as in the above case. 
However, a similar argument as in m Proposition 5.1] gives the following estimate with a 
log loss. 


W 3 (t, s)f(s)ds 


q. 


<h 1 \ log h | l/l L 2 . 

t,x' 


We complete the proof by gathering the above calculations to obtain 


Wi(t,s)f(s ) ds 


L l x , 


< G(h, a) 2 \f\\ L p> 

^ + r,J 


i = 1,2,3. 
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□ 


Proposition 3.2 (Case lb). Assume the notations in Case 1: 

IK hD Xl - a 1 {x,hD„j))u{h)\ vt < h\\u(h)\\ L2 , 

and 

II OAd - a 2 {x, aD x/ ))v{a)\ L2 < a \v(a) || l2 , 
where and <9|, f ,a 2 are positive definite. Suppose that Ch < a < h, where C is as i 


m 


Proposition \3.1\ Then 
where for n — 2, 


uv 


Lp 


F 2 , p (h) = 

and for n > 3, (n,p) (3, 2), 

and 


~ F n , P (h) \u \ L 2 \\v \\ L2 , 

h~for 2 < p < 3, 
h~ 1+ p for 3 < p < oo; 

F n , p (h) = h~ n+1+rL p- 


F^ 2 (h)—h 2 1 log h\ 2 . 


Remark. When cr « h, F n ^(h) « G ntP (h,a), and in this case Proposition 3.2 gives the L v 
bilinear quasimode estimates in full range. 

Proof. Note that since Ch < a < h the ratio cr/h ~ 1. Assume that ||m||l 2 = ||u||l 2 = 1- By 
scaling £ we obtain 

W(hD. X1 -d 2 (x,hD xl ))v \\ L 2 = 0(h), 

where 

a 2 (x,f) = ^a 2 (x, ^ . 

Therefore, we can treat both u and v as quasimodes at the same scale. We adopt the notation 
t — x i and write x G M as (t,x r ). Then consider the function 

f(t,x',z') = u(t,x')v(t,zf). 

We want to calculate the L p norm of / restricted to the submanifold x' — z'. Note that 

(h,D t - ai(t,x’,hD x f) - d 2 (t, z!,hD z ,))f(t,x’,z’) = 0 L 2(h). 

Therefore, we may directly apply the submanifold restriction estimates of m Theorem 1.7]. 
When n = 2, we are restricting from a 3-dimensional space to a 2-dimensional hypersurface, 
so we obtain 

h _2 + 2 t for 2 < p < 3, 

_1+ p for 3 < p < oo. 

For n > 3, we are restricting from a (2 n — l)-dimensional space to an n-dimensional subman¬ 
ifold, so we obtain 

F n , p (h) = h~ n+1+ P 

Note that as in in we must concede a log loss for the case (n,p) = (3,2), (in which case, we 
restrict a 5-dimensional space to a 3-dimensional submanifold). □ 


F 2 , p (h) = 
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3.2. Propagating in different directions (Case 2). Recall the notation x = (x\,x 2 ,x) 
and f = (6,6,0- 

Proposition 3.3. Assume the notations in Case 2: 

I (hD Xl - ai(x,hD X2 ,hDx))u(h)\\ L2 < h\u(h)\\ L2 , 

and 

II(cr D X2 - a 2 {x,aD xl ,aD s ))v(a )\\ L 2 < a \\v(a)\\ L2 , 
where <9^ and ^a 2 are non-degenerate. Then 

\\uv\\ LP 6 G U}P (h, cr) | u | | u | • 


Proof. Let 
and 


E h [u\ = ( hD Xl - ai(x, hD X2 ,hD x ))u(h), 


E a [v] = (cr D X2 - a 2 (x,aD xl ,aD x ))v(a). 

We set X\ = t\ and x 2 = t 2 and express u and v via the propagators 

( (hD t i - ai{s + h, x 2 , x, hD X2 , hD s ))U h (ti, s) = 0, 
\U h (0,s) = ld, 


and 


f ( aD t2 - a 2 (xi,s + t 2 ,x, hD Xl , hD x ))U a (t, s) = 0, 

{£4(0,6 = Id- 

Again we write £4(6) := £4(4,0) and £4(6) := £4(6,0) for notational convenience. Now 
we may write 


i 


u(ti,t 2 , x) = U h (ti)u(0, t 2 , x) + - £4(4 - Si, si)E h [u\ dsi, 


and 


So 


f*2 


v(t 1 ,t 2 ,x) = U a (t 2 )v(t 1 ,0 1 x) + - / U a (t 2 - s 2l s 2 )E (T [v\ ds 2 . 


a 


f*ti 


uv(t 1 ,t 2 ,x) = U h (t 1 )u(0,t 2 ,x)U a (t 2 )v(t 1 ,0,x) + -[£ ff (f 2 )ti(fi,0 1 i) / - s u s 1 )E h [u]ds 1 


ft 2 


H— U h (ti)u(0,t 2 ,x) / £4(6 - s 2 ,s 2 )E a [v\ ds 2 


a 


ah 


ftl rt2 


u h (t 1 - Si, si)E h [u\U a (t 2 - s 2 , s 2 )E a [v\ ds 1 ds 2 . 


'0 JO 


As in Proposition |3.1[ it is enough to prove 
£4(4)/(6, x)U a (t 2 )g(t i, x 


L ,. HG n , p (h,a)\f\ L , 


Using the parametrix construction, we have 

1 


U h (ti)f(t 2 ,x) = — (v&')bi(t 1 ,t 2 , X ,€ 2 ,€)f(s 2 ,y)ds 2 dyd£ 2 d£, 


and 


h n ~ 

1 


14 ( 66 ( 4,6 = / e AMt 2 M,x, vl ,v) s im <«'.’?»6 2 (^ lj t 2ja ; j7 y lj7 y)^( Sljt y)d Sl d u ,^ 1 d^ 


where 


<66(4, 6, z, 6,0 - ai(ti,t 2 ,x,d X2 <f> 1 ,'V s </>i) = 0, <6(0, t 2 , x, 6,0 = 66 + 6,0, 
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and 

- a 2 (t 2 , *i, a, §ei V*0 2 ) = 0, <6(0,6,^, 771,77) = + (.x,fj ). 

We will sometimes write t = (ti,£ 2 ) for notational convenience. Therefore, we need to study 
the A 2 (M 2 ) x L 2 (M 2 ( n ~ 2 )) —>• L P (M 2 ) x L p (M n_2 ) mapping properties of the operator Wh i<7 (t) 
given by 


(W h , a (t)h)(x) 


(ha] 


71—1 


h(si, s 2 , |/, w) d£id£ } dr] 2 dfjdsids 2 dydw. (3.22) 


To this end, we evaluate W/ ljCT (i)W^ (T (s): 

(W / / i „ (r (t)iy^ CT (s))/(x) = J W(t, x, s, z)f(s, z) dzds. 

Here, 


W(t, x , s, z) = 
where 

and 


(/nr) 2 !"- d 


r> s> £ C 2 , C, r/1,77, n, f) rfA, 
dA = drdydwd&dCd&dCdrjidfjdTidf, 


il’h.ait, r, s, x , 3/, w, z, 6, 6 C2, C, Vi, V, r i, T ) 

- (f>i(t 1 ,t 2 ,x,£ 2 ,£) -<f>i(s 1 ,s 2 ,z,C2,() ~r 2 (£ 2 - C2) - (y,€-() 


a 


h(h,t 1 ,x, 7 ] 1 ,ri) - 0 2 (s 2 ,si,z,ri,T) -r^ - n) - {y,y-r) 


We use stationary phase to calculate the (r 2 , y, (2, C, iq, in, Ti, t) integral. As the stationary 
point is always non-degenerate, we obtain 

W(t,x,s,z) = f e^ Mt i’ fe ^ 2 ’«)-^( s i’ S2 ’^ 2 ’d)6(t,s,x,f,6,x)de 2 ^ 


x —— / e i(Mt 2 M,x,vuv)-Ms 2 ,si,z,vimu tj 5> 5 771 , 77) drjidfj 

cr n J 

:= K h (t,s,x,z)K a (t,s,x,z) 

Now as in the proof of Proposition |3.1[ we analyse AT and K„ by studying the critical points 
of 


<t>i{ti,t 2 ,x,£ 2 ,€) ~ <i>i(si,s 2 ,z,£ 2 ,g) 

in (6,0 

(3.23) 

and 



<j) 2 (t 2 ,t 1 ,x,rii,fj) - 4 > 2 (s 2 ,si,z,rj 1 fj) 

in ( 771 , 77 ) 

(3.24) 

For both (3.23) and (3.24) to have critical points we must have both 



0 — (t 2 — s 2 ) + (x — y) (1 + 0(|s|)) + (A — s i)(<% 2 ,O ai ^(1^ — S D), 


and 


0 = (ti — si) + (x - y) (1 + 0(|s|)) + (t 2 - s 2 )(d {riun) a 2 + 0(\t - s|)). 

Since both |<9(£ 2j |\ai| and |9( r?1 ^)a 2 | are less than e, to have both critical points we require 

\t 2 - s 2 | < \ti — si|e, and |A - si| < |6 - s 2 |e. 
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Clearly, if e is chosen small enough this is impossible. Therefore, we are always able to 
integrate by parts in one of (^2,0 and (771,77). We split the kernel into two parts 

W(t,x,s,z) = Wi(t,x,s,z) + W 2 (t,x,s,z), 

where 

|^2 - S 2 | 


and 


W x (t,x,s,z) = x 


W 2 (t,x,s,z ) = 1 - x 


|*i - si 
1^2 — s 2 


W(t,x,s,z ), 


W ( t , x, s, z). 


,|*i - si | 

On the support of W\(t,x,s,z), we have \t 2 — S2 1 < |U — si| and therefore we cannot find a 
critical point in (771,77). So integrating by parts we obtain 


\K a (t,s,x,z)\ <a (n ( 1 + 


|*2 - S 2 | 


CT 


—N 


1 + 


\X — Z 


a 


-N 


and so using the trivial h b estimate on Kh we have 


|Wi 


(t,x,s,z)\<a-( n ~Vh-( n -V + + 


Using this estimate we obtain the L\ —> L 5 ° and L? —* L? norms of l-Uj (f, s) via Young’s 
inequality. Then we may use the Strichartz formalism as in Case 1 to get the L p ti —> 
norm of W\. 

Now on the support of W 2 we have 1 1 \ — s 1 < \t 2 — s 2 \ therefore there may only be critical 
points for K a but not Kh- This is the same as in the second regime in the proof of Proposition 
13.11 so we inherit those estimates here. □ 

4 . Sharpness of the L p bilinear estimates 

4 . 1 . Flat Model. We study the flat model, that is, localised quasimodes of the Laplacian 
in M n , to gain insight into sharp examples. 

Suppose that u is an L 2 normalised 0 12(h) quasimode of A K n. Then under Fourier trans¬ 
form 

where T\ t is the semiclassical Fourier transform defined as 

Then Th [ti] must be located near the sphere of radius 1 in the £-space. We create a family 
of quaismodes indexed by a which controls the degree of angular dispersion of £. Write 
£ = (r,oj) where u 0 G S n ~ l and set the coordinate system so that (u 0 corresponds with the 
unit vector in the £1 direction. Let 


X«(r,w) = 


1 if \r — 1| < /i, |u — ujq I < h Q , 
0 otherwise. 


Thpn Qpf 

&((,) = f h Jr,u) = r,u>). 

Note that is L 2 normalised. Now set 

U‘W = A- 1 (/t)W = (^75 
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is an L 2 normalised 0(h) quasimode of Ar™. We may write 


T*(x) = 


fo-l/2-a{n-l)/2-n/2 e ixi 

(2n) n / 2 


e UMSi-i)+(*,?)) Xa ^ d £' 


Note that if |xi| < eh 1 2,1 and \x'\ < eh 1 a for sufficiently small e > 0, the factor 

e £(®i(£i-l)+(z , ,£ , >) 

does not oscillate so in this region 

\T£(x)\ > ch~^ n - i)/ 2 +«("-i)/ 2 . 



x, 


Figure f. is localised so that is large in a h 1 2a x ( h 1 a ) n 1 tube 

Now suppose we have two semiclassical parameters a < h. For each we have a family of 
tubes T£ and Tf. We will construct sharp examples by studying the products of such tubes. 

4.1.1. L p bilinear quasimode estimates for large p. First we treat the case when p > 

We set a = fd = 0 and consider the product TfT^. 



sh 


-► 

*1 

Figure 2. To saturate the high p estimates we superimpose Tf and Tf 
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For |x| < £cr, we have 

\Tg{x)\ >a- (re - 1)/2 , and \T^x)\>h~^ 2 . 


So 


II T°,T { 


0 \\Lp ^ O' 


which therefore saturates the estimates of Theorem 


(n-i)/2^-( n -i)/2 a n/p _ C3 } l -(. n - 1 )/ 2 a - 

when p > 


1.3 


( n—l)/2—n/p 

5 

2(n+l) 
n—1 


4.1.2. L p bilinear quasimode estimates for small p in dimension n > 3 and midrange p in 
dimension two. We now construct saturating examples for 2 < p < in dimension 

> 3, these examples also saturate the estimates for 3 < p < 6 in dimension two. We choose 
f. Then we select the tube T° h such that 

a i-2a h = h ' 





Figure 3. To saturate the mid-range p estimates we superimpose Tfr and Tf h 
where ah is choose such that a 1 ~ 2o ‘ h = h 

Then for |a?i | < eh and \x'\ < ecr 1- ", we have 

\T 0 h (x)\ > h~ {n ~ 1)/2 , and \T° h (x)\ > a -(n-i)/ 2 +a (n- D /2 

So 

||Tq 1 t ct || > / 1 -( n - l )/ 2 cr -( n - l )/ 2 +“( n_l )/ 2 / 2 , l /p CT ( l - Q )( Ti - l )/p 

— (Jf l -3(n-l)/4+(n+l)/2p] l (n-l)/4-(n-l)/2p a .-(n-l)/2+(n-l)/p+a((n-l)/2-(n-l)/p) 

— ^</ 1 -3(n-l)/4+(n+l)/2p 0 .-(n-l)/4+(n-l)/2p 

This saturates the sharp L p bilinear quasimode estimates for 2 < p < in dimension 

n > 3 and 3 < p < 6 in dimension two. 
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4.1.3. IP bilinear quasimode estimates for 2 < p < 3 in dimension two. Here we set a — [3 — 
1/2 and consider the product Ti/ 2 (a)Ty 2 (h). 



Figure 4. To saturate the low p estimates we superimpose Ty 2 and TZ 2 


For \xf < e and \x 2 \ < ea 1 / 2 , we have 

\T°/ 2 0 )\> <M /4 , and |T 1/2 (ft)| >h- 1/t . 

So 

\\ T i/2 T i/2\\ LP > h~ 1/4 a~ 1/4 a 1/2p = Ch~ 1/4 a~ 1/4+1/2p , 
which saturates the estimates for 2 < p < 3 in dimension two. 


4.2. Spherical harmonics. In this subsection, we construct eigenfunctions on the sphere 


that saturate the L p bilinear eigenfunction estimates in Theorem 1.4 These examples also 


reflect the same behaviour as the flat model in the previous subsection. 

We first recall some standard facts about spherical harmonics as the eigenfunctions on the 
sphere. The spherical harmonics are the homogeneous harmonic polynomials restricted on 
the sphere; we use §13*. to denote the set of such functions with homogeneous degree k. For 
each u G §13*., u is an eigenfunction of — A §2 on § 2 with eigenvalue k{k + 1): 


—A§2 u = k(k + l)w. 

Therefore, the eigenfrequency A = \/k(k + 1) ~ k. The multiplicity of the eigenvalue k(k + l) 
is the dimension of its eigenspace, dim §33*, = 2k + 1. We use spherical coordinates 0 G [0,7r] 
and 6 G [0,27r) so that § 2 9 i= (sin 0 cos 6 , sin 0 sin 6 , cos 0 ). One can write the standard 
orthonormal basis of §33*, as 


= C„, m PT(c os^e*"*, 

in which Ck,m is the L 2 normalisation factor, and P™ is the associated Legendre polyno¬ 
mial. We remark that the phase of Yjf is constant on each longitude (i.e. 9 = const), and 
their modulus is constant on each latitude (i.e. 0 = const). Two special cases of spherical 
harmonics are as follows. 


(i). m = 0: Z k = Y, 


o 


points 0 = 0 and 0 
(1.1) for large p: 


are called zonal harmonics. Z k concentrates on the two antipodal 
= n. They achieve the maximal norm growth in Sogge’s L p estimate 

i-ZfcllLp ~ / c 2 C 1 / 2 - 1 /p)- 1 / 2 , 6 < p < oo. 


If 101 ^ k x , i.e. within a k 1 neighborhood of the north pole, |Z*,| > k 1 ^ 2 . Similar 
estimates holds also around the south pole. 
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(ii). m = ±k: Q±k = Y± k are called highest weight spherical harmonics or Gaussian beams. 
Q±k concentrate in a fc _1//2 neighborhood of the equator f = tt/2 and achieve the 


maximal norm growth in Sogge’s L p estimate (1.1) for small p: 


\\Q±k\\Lp « k {1/2 ~ 1/p)/2 , 2 < p < 6. 

Notice that Q±k decreases exponentially away from the concentration and ||Q±fc||L° 


To construct the sharp examples for L p bilinear eigenfunction estimates in Theorem 1.4 


we 


divide the range of p into large p {p > 6), midrange p (3 < p < 6), and small p (2 < p < 3). 

4.2.1. L p bilinear eigenfunction estimates for p > 6. Let |_AJ be the largest integer that is 
smaller than A. Write u = Zy \j and v = Z . Then 

\u{f,0)\> if |0|<eA _1 , 

and 


\v{<f>,9)\ > P*, if \(f\<£p l . 


So 


\uv\\ LP > A 2 /i 


1 _ 2 
2 P. 


which saturates the estimate in Theorem 1.4 when p > 6. 


4.2.2. L p bilinear eigenfunction estimates for 2 < p < 3. Let u = <3|aj and v — Q]ji\- Th en 

7r 


and 


So 


K<M)I>at if 

if 


< eA" 


< £jl 


\\ uv \\lp > 

which saturates the estimate in Theorem 1.4| when 2 < p < 3. We remark that on § n when 
n > 3, Burq, Gerard, and Tzvetkov (BGT4j proved that two zonal harmonics saturate the 
L 2 bilinear eigenfunction estimates in Theorem 


1.1 


(modulo log when n — 3). 


4.2.3. L p bilinear eigenfunction estimates for 3 < p < 6. For the smaller eigenfrequency A, 
we let u — Zyx j- Then 

Hf,0)\>XZ if |0| < sA -1 . 

Let a G [0,1/2). We set 

! 1 = \t=2£, (4.1) 

and construct the eigenfunction v such that 




1 —o: 
2 


if |aq| < e\ 1 = ep i 1 20 ^ and \x 2 \ < ep (1 a \ (4.2) 

3 (xi, x 2 ,x 3 ) = (sin (p cos 9, sin0sin 9, cos 0). One sees imme- 


Recall that in our notation, 
diately that this is the eigenfunction version of the flat modal in §|4.1.2| In view of (|4.lh 


I II \ \ i 1=2 , -1 -1=2 

\uv\\ L p £ A 2 p 2 A pp p 


■k_3_ 1_JL 
> A 4 2 p p 4 2 P. 


To construct the required spherical harmonic v, we use linear combination of Gaussian beams. 

Set k — [p\. According to the different propagating directions of Qk and Q-k, we say that 
Qk has the north pole <f = 0 as its pole and Q_k has the south pole </> = tt as its pole. Since 
A is rotational invariant, given any point p e § 2 , we can find a Gaussian beam with pole p 
by rotating Q k - Two Gaussian beams with the same pole only differ by a phase shift. (See 
[ Hal, Lemma 8] for more details.) Hence, a Gaussian beam in SH fc is uniquely determined by 
its pole and the phase of some point on the sphere. 
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Next we need a theorem in [Ha], which shows that for a family of well-separated poles 
{Pj}]Li C § 2 , the corresponding Gaussian beams {qj}Jk 1 C SHI*, are almost orthogonal. This 
enables us to estimate the L 2 norm of a superposition of well-separated Gaussian beams. 


Lemma 4.1. There exists a positive number d > 0 such that the following statement is true. 
For any family of well-separated poles {pj}jLi C § 2 satisfying 

dist(pj,pj) > d ■ k ~2 for all i,j, 

the eigenvalues of the Hermitian matrix (( q^qj)) are all in [1/2,2], Here, qj is the Gaussian 
beam with pole pj and (•, ■) denotes the inner product in L 2 (§ 2 ). 

The proof can be found in m Section 3], 

Let m = [k 1 /' 2 ~ a \. We choose the poles {pj}fLi around 6 = 7t/2 on the equator 0 = 7t/2 
with distance between each pair satisfying dist (qi,qj) = d ■ fc -1 / 2 . Then all pj fall into a 
0o-neighborhood of 6 — 7t/ 2 on the equator, where 

6q = m ■ d ■ k~ 2 ~ k~ a . 


Without loss of generality, assume that p\ = (0,1,0), that is, 0 = 7 t/ 2 and 6 = 7t/ 2 in 
spherical coordinates. Then the corresponding Gaussian beam qi concentrates around the 
great circle defined by the equation 9 = 0,7r. Any other corresponding Gaussian beam qj 
concentrates on the great circle which has an angle < 9 0 with the one of q \. And q'-s intersect 
at the north pole p n . Write 

m 

W =—^ qj . 


m 


3 =1 


Then thanks to Lemma 4.1, we have 


Mil* = — 

m 


i=i 


1 

m 




9.7 ^ 


< 2. 


L 2 


2=1 


3=1 


We now set the phase of all q'-s at the north pole to be e lk0 = 1. Hence, 

1 


w(p n ) = ~/=y^dj(Pn) >k 2 
\ m 


i=i 


We need to show that the above lower bound holds in the neighbourhood of p n 
S = {x = (xi,X 2 ,x 3 ) G § 2 : |xi| < £k~ ( - 1 ~ 2a ' > and \x 2 \ < £k~ ( ' 1 ~ ol ' > } 


(4.3) 


for sufficiently small £. Notice that S falls into the concentration tube of every q 3 since 
a E [0,1/2) and \x 2 \ < e/c _( ' 1_Q ^ < k _1//2 . Therefore, £x x e S, then 

\Qj( x )\ ~ f° r a U .?• (4.4) 

To determine the phases of qj(x) for j = 1,..., m, we let pj be the new north pole and denote 
dj(x) as the longitudinal difference of x and p n in this new coordinate system. Hence, the 
phase difference of q 3 (x) and qj(p n ) is e lkdj ^ x \ therefore the phase of q 3 (x) is e tkdj ^ since the 
phase of qj(p n ) is set to be 1. 

We observe that for any x G S, 

\di(x) — dj(x )| < ■ 9 0 < £k~ l . 

Thus, by choosing e small, the real part 

3? > 1 for all ? = i, m . 
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Hence, 


1 1 

w(x) = — | qj(x))\p ikd ^ x) = - * ikdl{x) 

j =i 


m 


m 




D ik(dj(x)—di(x)) 


and then in view of (|4.4[) 


j =i 


3? (e- ifcdl(a:) w(a;)) = -= |g,(x)|3? (e ifc ( dj '( x )- d i( x ))) > k 

\/ TXl 


l-OL 

2 


i=i 


Let v = w/\\w\\ L 2. Then by (4.3), 

1 — ck 

|n(x)| > k for all x G S, 


which gives the required spherical harmonic in (4.2). 


5. Appendix: Semiclassical analysis 


In this appendix, we provide the background on semiclassical analysis that is used in this 
paper. We refer to Zworski |Zwj for a complete treatment in this subject. 

Semiclassical analysis provides a valuable framework in which to consider high frequency 
problems. The key idea is to scale a fixed frequency A —>• 1. This scaling gives rise to a 
semiclassical Fourier transform (with X = 1/h) 

= (27r/ ^ n/2 J e~^u{x)dx, 

and its inverse 

Note that as with the standard Fourier transform there is some variance in the literature as 
to the pre-factor. We choose (27T h)~ n ^ 2 as this is exactly the right factor to preserve the L 2 
norms. That is, 

\\FhU\\ L 2 = HI £2 • 

Under this scaling the standard relationships between differentiation and multiplication are 
retained (although they now also include a scale factor). That is, 

F h [hD Xi u] = iiT h [u\ and F^ l [xif] = -hD ii J r h ; 1 [f}. 


Therefore constant coefficient differential operators can be written as 

1 


a(hD)u(x ) = 


(27rh) T 


e' h{x v ’°a(£)u(y)d£dy. 


Naturally we want to extend to semiclassical pseudodifferential operators a(x , hD)\ 

1 


a(x, hD)u(x) = 


(2tt h) r 


e l(x y& a ( Xi £)u(y)d£dy. 


Therefore we must define some symbol classes for the symbols a(x,£) to live in. 

Like the standard pseudodifferential calculus our symbols will be classical observables on 
phase space. Suppose that (M, g) is a Riemannian manifold. An element in the cotangent 
bundle T *M is denoted as (a;, £) with x G M and £ G T* M. We write |£| x as the induced 
norm of £ G T*M by the Riemannian structure g. When we work locally, as we often do, we 
may associated T *M with patches of M 2n . In this paper we concern ourselves with functions 
that are semiclassically localised. 


Definition 5.1 (Semiclassical localisation). Let K be a compact subset of T* M. We say 
that a tempered family {u(h)} is localized to K in the phase space if there exists a function 
X G C^(K) for which 

u(h) = x( x , hD)u(h) + Os(h°°). 
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This means that we may assume that any symbols we use have compact support. In 
the standard pseudodifferential calculus, symbols with compact support are smoothing and 
therefore do not give rise to interesting mathematics. However in the semiclassical calculus 
the h scaling means that even smooth and compactly supported symbols can indeed give rise 
to very interesting mathematics. 

We need to define a sensible class of symbols and a procedure to quantise them. That is, 
to associate them with a semiclassical pseudodifferential operator. 


5.1. Symbol classes. In the analogy to the standard calculus one may define symbols a(x, £) 
that on compact subsets K of T *M 

sup |S£dfa(a;,£)| < C a: p, K (£) m_l/31 

x£K,£eT£ M 

for some C a> p t K independent of h. However these classes don’t capture the behaviour of the 
symbol with respect to the semiclassical parameter. To motivate our choice of symbol class, 
consider the standard symbols 

a(z,f;A) = (l + l£D T p(y) • 

where p G Co°(M) has support near 1. Hence, these symbols are localised near |£| = A. We 
have that 

a(x, D)u(x) = -i- J (l + |S| 2 ) ¥ p f ■ u(y) d^dy, 

re-scaling £ —>■ A£ and setting A -1 — h we have 

a(x, D)u(x) = t—— [ e^ x ~ y ’® 

y y ’ (2? rh) n J 

where 

a h (x,&h) = ^1 + p( |£|). 

Note that the symbol a,h(x, £) is localised near |£| = 1 so 


£ 2 

1 + — 
h? 


p(\t,\) u (y) d&y = a h(x, hD)u(x), 


\a h (x,^h)\<h- m . 


Therefore the level of smoothness (captured by the order m) of a symbol in the standard 
calculus translates to decay in powers of h when we consider it semiclassically. Therefore we 
define a set of symbol spaces of compactly supported symbols to reflect this relationship. 


Definition 5.2. Let m G M. We define the symbol classes 

S m (M) = {a G C™(T*M) | \8 a a\ < C a h~ m }. 

(1) If m — 0, we denote S m ( M) by S'(M). 

(2) We denote *S-°°(M) = n meR S m (M) and = U meR S m (M). 

If a(x, £; h ) is independent of h we write it as a(a;, £). We refer the reader the Zworski [ZwJ 
for definitions of symbols that are not compactly supported or that lack of smoothness as 
/i -> 0. 


5.2. Semiclassical pseudodifferential operators. Having an appropriate class of sym¬ 
bols we can now associate every symbol with a semiclassical pseudodifferential operator (i.e. 
TDO). As in the standard calculus there is some choice to the quantisation procedure. 


Definition 5.3 (Standard and Weyl quantisations). Given a G S m (R n ), we define 
(i). the left (or standard) quantisation as 


a(x, hD)u(x) 


1 

(27t h) n 


e h( x a{x,£\ti)u{y) dt^dy for u G 5(M n ); 
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(ii). the Weyl quantisation as 
a w (x, hD)u(x) = 


(2nh) r 


Mx-y£) 


a 


'R 2 n 


x + y 


, £; h ) u(y) d^dy for u G d>(R n ). 


If the symbol is instead in S' m (M) for (M, g ) a smooth Riemann manifold we may keep 
this definition replacing M 2n with T*M. 

While there are other choices of quantisation these two choices are often the most useful. 
The Weyl quantisaton has the nice property that a w (x, hD) is self-adjoint if a is a real-valued 
symbol. The standard quantisation has the nice property that 

a(x, hD)u = (a(x, £; h)F h [u \). 

Remark. If a(x,£\ h) G 5 m (M) then one can show (by almost orthogonality) that 

\\a(x,h,D)\\ L ^ L 2 <h- m . 

In particular, if a(x,£;h) G S'(M) then the L 2 —>• L 2 mapping norm of a(x,hD ) is bounded 
independent of h. 

We now define the set T m (M) of semiclassical pseudodifferential operators with symbols 
in and establish the correspondence of A G T 7n (M) and its semiclassical principal 

symbol a. The correspondence is one-to-one modulo lower order terms. Denote 

a = a(A) : T m (M) —>■ S m (M)/ S m ~ l (M), 

and its right inverse, a non-canonical quantisation map for a G S' m (M): 

A = Op h (a) : S m ( M) ->■ T m (M). 

&(A) is called the principal symbol of A. It is modulo S' m-1 (M) unique under change of 


quantisations and change of local coordinates. In the same fashion in ( 5.1 

(1) if m = 0, we denote T m (M) by T(M); 

(2) we denote 'h _0 °(M) = and ^°°(M) = U mS R\I /m (M). In this context the 

elements in T~ 0 °(M) are referred as smoothing operators; 

We sometimes abuse notation somewhat and neglect to show the cutoff functions that 
provide the compact support. For example we may work with the Laplacian and refer to its 
symbol as |£| 2 when in fact in our localised setting we are really working with the operator 
p(x, hD) with symbol 

P(x,0 =p(|£U)l£lx 

where p G C*^°(M). 

The usual operations involving semiclassical pseudodifferential operators are as follows. 
Let A G T m (M) and B G T m, (M). 

(1) Let A* be the adjoint operator of A in L 2 (M). Then 


( 2 ) 


( 3 ) 


a(A*) = a (A) +r(x,£;h), 
where the error term r(x, £; h) G 5 ,,7l_1 (M). 

cr(AB) = a(A)a(B) + r{x, & h), 
where the error term r(x, h ) G S m+m ^ 1 (M). 


(5.1) 


(5.2) 


^([A B \) = -ih{a(A), a{B)} + hr(x, h), (5.3) 

where the error term r(x,C,] h ) G S m+m _1 (M) and {•, •} denotes the Poisson bracket 
defined by 

, , da db da db 

dx dl C dC, dx 
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Note that like the standard calculus, error terms from commutation and composition are 
of lower order. In the semiclassical setting this means that each of these error terms is one 
factor of h better than the first term. 

We now briefly discuss some standard properties and estimates from semiclassical analysis 
relevant to this paper. 


5.3. Elliptic operators. We say that a G 5(M") is elliptic if |a(a;,£)| > c > 0 for all 
(x,£) G T*R n . If a is elliptic, then there exists b G S(R n ) such that b(x, hD ) is the inverse of 
a(x, hD) modulo a smoothing operator. In fact, this assertion holds in the microlocal sense: 


Theorem 5.4 (Inverses of elliptic TDOs). Let x € 5(R”). If a G S^R* 1 ) satisfies |a(:r, £)| > 
c > 0 for all (x,£) G suppy. Then there exists b G S( R n ) such that 

b(x, hD)a(x , hD)x(x, hD) = x( x , hD) + 0 L 2^ L 2(h°°), 

and 

a(x, hD)b(x, hD)x(x, hD) = x(x, hD) + O i 2 _^ Z/ 2 (/i°°). 

These equations also hold if we use Weyl quantisation instead of the standard one. Moreover, 
if x £ C£°(T* R n ) ; then we can replace 0 L 2^. L 2(h°°) by Os'^s(h°°). 


The proof of Theorem 5 A is similar to the analogous statement for the standard calculus. 
However, in this setting we express everything in powers of h. If a(x,£ m , h) is bounded away 
from zero, then we may define bo(x, hD) by its symbol b(x,£] h) = l/a(x, £; h). Since a(x, £; h) 
is bounded away from zero bo(x,£; h) G S^R" - ) = S°(M. n ) Then the composition formula (5.2) 
gives that 

a(b 0 (x, f; h)a(x, f; h)) = 1 + r x (x, f; h), 

where ri(x,£]h) G S' m ~ 1 (R n ), that is, it is one order of h better than a(x,£]h). Then we 
may set 

7 , c M ri(x,f;h) 
b ^ h )= - r. h) ’ 

and the composition formula ( 5.2[ ) gives 

(b 0 (x, hD) + hbi(x, hD)) a(x, hD) = 1 + r 2 (x, hD), 

where r 2 (x,f > ] h) G S m ~ 2 (W l ). We can continue this process to produce 

N 


b N (x, f; h) = ^2 h l bi(x, f; h) 


i =0 


such that 


b N (x, hD)a(x, hD) = 1 + rjy(x, hD) 

with rjsr(x,£',h) G S m ~ N . Therefore modulo an 0(h°°) term we may produce an inverse for 
a(x,£;h). 

The localisation property alone is enough to prove a range of L p estimates (see e.g. [KTZl 
Lemma 2.2]) 


Theorem 5.5 (Semiclassical L p estimates). If a G S'(R), then 

a(x, hD) = O : L p (R n ) —> L q ( R n ), 

in which 1 < p < q < oo. 


An immediate consequence of Theorem 5.5 is that if u = u{h) is a family of localised 
functions then (as we may locally teat T*M as R 2n ) 
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5.4. Quasimodes. Let A(h) E T m (M). We define quasimodes as approximate solutions to 
A(h)u(h ) = 0. An 0 L 2 (h a ) quasimode is a function u E L 2 such that 

\\A(h)u(h )\\ L 2 < h a \\u(h)\\ L 2 . 


The finest quasimode resolvable by semiclasscical analysis is 0^2 (h°°) that is 


\\A(h)u(h)\\ 


L 2 


<h N 


u\ 


L 2 


for any N. 

We are particularly interested in the case when a = 1, that is, (h) quasimodes. These 
quasimodes behave well under localisation. Suppose x( x >0 is a smooth compactly supported 
function then x{ x i£) £ S'(M) (we should think of this function as having small support in 
phase space). By the composition formula (5.2), 

A(h)x{x, hD)u = a(x, hD)x(x, hD)u = xi x i hD)a(x, hD)u + r(x, hD)u, 


where if a(x,£; h ) G S^M), then r(x,£; h) E S m 1 (M). If a(x,£; h) E 5(M), which is often 
the case, then r(x,£] h) E S'^ 1 (M) and 


\\r( x ,hD)\\ L2 ^ L2 < h. 

Therefore if u is an 0^2 (/i) quasimode of A(h), the function v = x( x , hD)u is also an 0^,2 (/i) 
quasimode of A(h). This means that we may always decompose phase space into small size 
patches and work on each patch separately. Finer quasimode order is not necessary preserved 
under localisation making the treatment of such function technically difficult. Finally when 
u is an 012(h) quasimode and a(x, £; h) E S'(M) we may ignore the dependence of the symbol 
on h. We write 

a(x,^\ h ) = a o(x, £) + lm(x , £; h), 

where ao(x, £) = a(x,£, 0) and a E S. Therefore if u is an 0^2 (h) quasimode of a(x,hD ), 
then it is also an Oi?(h) quasimode of a 0 (x, h,D). For this reason we only study quasimodes 
of operators associated to symbols a(x,£) independent of h. 


Example. Let P(h ) = — h 2 A — 1. Consider 

A Uj = A 2 jUj, where |wj|| i2 = 1- 

This family of eigenfunctions {u(hj)} with hj = XJ 1 is an quasimode of any order (i.e. 
0 L 2(h°°) quasimodes). 

Now define the spectral clusters as linear combinations of eigenfunctions in a spectral 
window with fixed length: 

u v = 'y ^ c j u j- 

Aj G[/i,/x+l] 

Then u(h) = ui/h is an 0^2 quasimode to P(h). The same conclusion also holds if one 
considers the smoothed spectral cluster (i.e. Sogge operator) 

XfiU = x [V -A - pju with x e C“(M). 

That is, u(h) = Xi/h u is a family of O l 2 quasimodes of P(h). See | Zwl §7.4.1] for more 
discussions on quasimodes. 


It is often instructive to look at the singular structure of families of functions. This is 
usually done through the semi classical wavefront set. The semiclassical wavefront set WF^(u) 
of a tempered family u(h) is the complement of the set of points (x, £) G T*M such that there 
exists a E C£°(T* M) with support sufficiently close to (x, £) with 

a(x, hD)u(h ) = Os(h°°). 
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Let A £ T m (M) with principal symbol a £ S m ( M) and u(h) be a Os(h°° ) family of quasi¬ 
modes, i.e. 

A(h)u(h) = O s (h°°). 


Then 


WF h (u) C {(x,0 G T* M : a(x,Q = 0}. 


Example. The eigenfunctions (as a tempered family) {uj}fL 0 = {u(/ij)}jL 0 are the solutions 
to 

P(h)u(h) = 0, 

where P(h) = —h 2 A — 1 has symbol p(x,£) = — 1. Then the semiclassical wavefront set 

WF h (u) is contained in p _1 (0) = {(:r,£) £ T*M : |£| x = 1} = S*M, the cosphere bundle. 


As stated above we intend to work with localised symbols. To translate our results to 
eigenfunctions we need to know that they are indeed semiclassically localised. 


Example. The eigenfunctions {uj}fL 0 = {u(hj)}fL 0 admit localisation property: Note that 
S* M is compact, one only needs to choose y £ Co°(M) such that supp x S* M and equals 
1 around S'*M, then 

(1 - x(x, hjD))u(hj) = O s (hf), 
since 1 — y G Qj°(M) and supp (1 — y) D WFa(u) = 0. 


5.5. Evolution equation and semiclassical Fourier integral operators. Let P £ T m (M) 
with principal symbol p. A key intuition in the study of quasimodes of P is to link their 
behaviour to the classical flow given by 

f x(t) = 

\£(t) = -d x p(x,£). 

We often think of a quasimode as being comprised of small wave-packets tracking along trajec¬ 
tories of the classical flow. To make use of this intuition we need to understand semiclassical 
evolution equation. 

Consider the inhomogeneous semiclassical evolution equation. 

J ( hD t + A(t))u(t, x) = f{t, x), for (t, x) £ M x M n , 

1 w(0) = u 0 . 


Here, A(t) G S(Ix W 1 ) is a family of TDOs with t as the parameter. The principal symbols 
a(t,x, £) = a(A(t)) are real-valued and independent of h. We may use Duhammel’s principle 
to reproduce u in terms of the propagator U(t) and the inhomogeneity. We represent U(t) 
by an h°°- approximate propagator, that is, a solution to 


(hD t + A(t))U(t)u = 0 L 2(h°°), for |f| < T, 
U(0)u = u. 


(5.4) 


We find a parametrix solution to (5.4) where U(t) is a semiclassical Fourier integral operator. 
Indeed for T small, 

1 f 


U(t)u(x) = 


~$ri ./TRm 


(27 ih) 

where E(t) = 0(h°°) : S' —> S, ip satisfies 


eh W>(h®,£) x ^- h)u(y) d£dy + E(t)u(x), 


d t (j)(t, x, rj) + a{t, x, d x (j)(t, x, £)) = 0, 

<f>(0,x,C) = (x,0, 


and 


b(t,x,i] h) £ C£°(M x T*R n x M). 
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Sketch of proof. We seek a parametrix solution of the form 

1 


U(t)u(x ) = 


(2nh) 7 


e fcW>(Lz,0 h)u(y ) d£dy, 


where 

0(0, x, 0 = (x, f) and 6(0, x, f; h) = 1. 

Note that this dehnitely satishes the initial conditions. Then 

hD t [U{t)u(x)\ = ^ h y n J [Mt, x, O b (t, x, & h) + hb t (t, x, f)]u(y) d£dy, 

and using the method of stationary phase 

1 


a(x,hD)[U(t)u(x)] = 


(2n h) 7 


Therefore, if 
then 

To improve the error write 


{y£}) [a(t,x,V x (f))b(t,x,£; h)+hr(t,x,£,; h)]u(y) dfdy. 
0t(C x, f) + a(t, x, V x .0) = 0, 


\\(hD t + A(t))U(t)\\ l 2 ^l 2 < h. 


N 


b(t, x, £; h) = 1 + '^2b j (t,x,£\h). 

3 = 1 

Then by solving transport equations for each of the bj we may improve the error to h N for 
any N. Full details of the proof can be found in jZwl Chapter 10]. □ 
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